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Intrinsic Localized Modes (ILM) (or Discrete Breathers, DB) are localized oscillatory modes 
known to occur in atomic or molecular chains characterized by coupling and/or on-site potential 
nonlincarity. Quasi-crystals of charged mesoscopic dust grains (dust lattices), which have been 
observed since hardly a decade ago, are an exciting paradigm of such a nonlinear chain. In gas- 
discharge experiments, these crystals are subject to forces due to an externally imposed electric 
and/or magnetic field(s), which balance(s) gravity at the levitated equilibrium position, as well 
as to electrostatic inter-grain interaction forces. Despite the profound role of nonlinearity, which 
may be due to inter-grain coupling, mode-coupling and to the sheath environment, the elucidation 
of the nonlinear mechanisms governing dust crystals is still in a preliminary stage. This study is 
devoted to an investigation, from very first principles, of the existence of discrete localized modes 
in dust layers. Relying on a set of evolution equation for transverse charged grain displacements, 
we examine the conditions for the existence and sustainance of discrete localized modes and discuss 
the dependence of their characteristics on intrinsic plasma parameters. In addition, the possibility 
of DB stabilisation via an external force is discussed. 

PACS numbers: 52.27.Lw, 52.35.Fp, 52.25.Vy 
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I. INTRODUCTION 

A variety of linear and nonlinear collective effects are known to occur in a dust-contaminated plasma (dusty plasma, 
DP) Q and relative theoretical research has received new impulse, since roughly a decade ago, thanks to laboratory 
and space dusty plasma observations. An issue of particular importance in DP research is the formation of strongly 
coupled DP crystals by highlycharged dust grains, typically in the sheath region above a horizontal negatively biased 
electrode in experiments 0, 0] ■ Low- frequency oscillations are known to occur Q in these mesoscopic dust grain 
quasi-lattices in the longitudinal (~ x, in-plane, acoustic mode), horizontal transverse (~ y, in-plane, shear mode) 
and vertical transverse (~ z, off-plane, optic-like mode) directions. 

Various types of localized (nonlinear) excitations are known from solid state physics to exist in periodic chains 
(lattices)oi interacting particles, in addition to propagating vibrations (phonons), due to a mutual balance between 
the intrinsic nonlincarity of the medium and dispersion. Such structures, usually investigated in a continuum approx- 
imation (i.e. assuming that the typical spatial variation scale far exceeds the typical lattice scale, e.g. the lattice 
constant 7*o), include non-topological solitons (pulses), kinks (i.e. shocks or dislocations) and localized modulated 
envelope structures (envelope solitons), andgeneric nonlinear theories have been developed in order to investigate 
their relevance in different physical contexts j3j- In addition to these (continuum) theories, which deliberately sacrifice 
fS |. discreteness in the altar of analytical tractability, attention has been paid since more than a decade ago to highly 
localized (either stationary or propagating) vibrating structures [e.g. discrete breathers (DBs), also widely referred 
• *-H . to as intrinsic localized modes (ILMs)], which owe their very existence to the lattice discreteness itself. Thanks to a 
few pioneering works 0, |H Q>1S an d a number of studies which followed, many aspects involved in the spontaneous 
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formation, mobility and interaction of DBs are now elucidated, both theoretically and experimentally; see in Refs. 
|9l lTol,ITH | for a review (also see Refs. 0,0, with reference to this study). 

Despite the fact that nonlinearity is an inherent feature of the dust crystal dynamics (either due to inter-grain 
electrostatic interactions, to the sheath environment, which is intrinsically anharmonic, or to coupling between different 
degrees of freedom), our knowledge of nonlinear mechanisms related to dust lattice modes still appears to be in a 
rather preliminary stage today. Small amplitude localized longitudinal excitations (described by a Boussinesq equation 
for the longitudinal grain displacement u, or a Korteweg-deVries equation for the density du/dx) were considered in 
Refs. |l lj and generalized in Ref. |15| . The nonlinear amplitude modulation of longitudinal and transverse (vertical, 
off-plane) dust lattice waves was recently considered in Refs. and 0,0 (also see 0), respectively. In 

fact, all of these studies rely on a continuum description of the dust lattice. On the other hand, the effect of the 
high discreteness of dust crystals, clearly suggested by experiments 0,0,0,0], may play an important role in 
mechanisms like energy localization, storage and propagation and thus modify the crystal's dynamical response to 
external excitations (in view of DP application design, e.g.). To the very best of our knowledge, no study has been 
carried out, from first principles, of the relevance of DB excitations with respect to dust lattice waves, apart from 
a preliminary investigation (restricted to single-mode transverse dust-breathers), which was recently presented |25| . 
This text aims in making a first analytical step towards filling this gap, by raising a number of questions which have 
not been addressed before. This study is neither exhaustive nor complete; it will be complemented by forthcoming 
work. 



II. THE MODEL 



We shall consider the vertical (off-plane, ~ z) grain displacement in a dust crystal (assumed quasi-one-dimensional: 
identical grains of charge q and mass M are situated at x n = n r*o, where n = ...,—1,0, 1,2,...), by taking into account 
the intrinsic nonlinearity of the sheath electric (and/or magnetic) potential. The in-plane (longitudinal, acoustic, 
~ x and shear, optical, ~ y) degrees of freedom are assumed suppressed; this situation is indeed today realized in 
appropriate experiments, where an electric potential (via a thin wire) 0] or a coherent light (laser) impulse 0,0,0 
is used to trigger transverse dust grain oscillations, while (a) confinement potential(s) ensure(s) the chain's in-plane 
stability. 



A. Equation of motion 

The vertical grain displacement obeys an equation in the form [l8l Il9j 

d 2 S z dS z 

, 2 " + v —^~ + u l ( Sz n+i + dzn-i - 2 8z n ) + Sz n + a (Sz n ) 2 + (3 (Sz n ) 3 = , (1) 
where Sz n (t) = z n {t) — zq denotes the small displacement of the n— th grain around the (levitated) equilibrium position 

1/2 

Zq, in the transverse (z—) direction. The characteristic frequency uiq = [— q& (r$) / {Mtq)\ results from the dust 
grain (electrostatic) interaction potential 3>(r), e.g. for a Debye-Hiickel potential 0|: $_o(r) = (q/r) e~ r ^ D , one 



has: D = q 2 /{Mr^) (1 + r /\£,) exp(— r /A£>) , where Ad denotes the effective DP Debye radius ll|. The damping 
coefficient v accounts for dissipation due to collisions between dust grains and neutral atoms. The gap frequency uj g 
and the nonlinearity coefficients a, [3 are defined via the overall vertical force: 

F(z) = F e/m - Mg a -M[u 2 g Sz n + a (Sz n ) 2 + (3 {Sz n ) 3 } + 0[(Sz n ) 4 ] , (2) 

which has been expanded around zq by formally taking into account the (anharmonicity of the) local form of the 
sheath electric (follow exactly the definitions in Ref. [Tgj n °t reproduced here) and/or magnetic 0] field(s), as well 
as, possibly, grain charge variation due to charging processes Recall that the electric/magnetic levitating force(s) 
F e / m balance(s) gravity at zq. Notice the difference in structure from the usual nonlinear Klein-Gordon equation used 
to describe one-dimensional oscillator chains — cf. e.g. Eq. (1) in Ref. 0: TDLWs ('phonons') in this chain are 
stable only in the presence of the field force F e / m . 

For convenience, we may re-scale the time and vertical displacement variables over appropriate quantities, i.e. 
the characteristic (single grain) oscillation period uj^ 1 and the lattice constant ro, respectively, viz. t = u~ x t and 
Sz n = roq n ; Eq. 0]) is thus expressed as: 

d 2 q 

+ e{q n+1 + q n - 1 -2q n )+ q n + a' q 2 n + f3' & = , (3) 
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where the (dimensionless) damping term, now expressed as (v /u) g )dq n / dr = v'q ni will be henceforth omitted in the 
left-hand side. The coupling parameter e = ujq/lu 2 measures the strength of the inter-grain interactions (with respect 
to the single-grain vertical vibrations); this is typically a small parameter, in real experiments (see below). The 
nonlinearity coefficients are now: a' — ar$/u 2 and [3' = ^r^jus^. 

Eq. (0 will be the basis of the analysis that will follow. Note that the primes in a! and (3 1 will henceforth be 
omitted. 



B. The model Hamiltonian 



In order to relate our physical problem to existing generic models from solid state physics, it is appropriate to 
consider the equation of motion £Q) as it may be derived from a Hamiltonian function, which here reads: 



H = 



N 

E 



£r + "<«> 



j 



(4) 



Here, pj obviously denotes the (classical) momentum pj — M<jj. 
environment, is of the form: 



The substrate potential, related to the sheath plasma 



(5) 



The coupling parameter e takes small numerical values (cf. below), accounting for the high lattice discreteness 
anticipated in this study. The minus sign preceding it denotes the inverse dispersive character of (linear excitations 
propagating in) the system; see the discussion below. Upon setting e — > — e, the 'traditional' (discretized) nonlinear 
Klein- Gordon model is recovered p8j| . 

It should be noted that both experimental [2lJ and ab initio (numerical) [2^ studies suggest that dust crystals 
are embedded in an nonlinear on-site (sheath) potential V, in the vertical direction, which is (possibly strongly) 
asymmetric around the origin, i.e. not an even function of qj. This implies a finite value of the cubic anharmonicity 
parameter a, thus invalidating models involving even potential forms - e.g. V(qj) ~ ?f/2 + /?9^/4 - in our case. 



III. LINEAR WAVES 



Retaining only the linear contribution and considering oscillations of the type, Sz n ~ exp[i (knro — tot)] + c.c. 
(complex conjuguate) in Eq. Q, one obtains the well known transverse dust lattice (TDL) wave optical-mode-like 
dispersion relation 

c 2 =^-4c 2 sin 2 ^), (6) 

i.e. 

Cj 2 = 1 - 4e sin 2 (fc/2) . (7) 

See that the wave frequency u) = oju> g decreases with increasing wavenumber k = 2ir/X = k/r (or decreasing 
wavelength A), implying that transverse vibrations propagate as a backward wave: the group velocity v g = u>'(k) and 
the phase velocity u> p h = tu/k have opposite directions (this behaviour has been observed in recent experiments). The 
modulational stability profile of these linear waves (depending on the plasma parameters) was investigated in Refs. 
|l8l Il9j|. Notice the natural gap frequency uo(k = 0) = uj g — uj m ax, corresponding to an overall motion of the chain's 
center of mass, as well as the cutoff frequency Lo m in — (w 2 — 4WQ) 1 / 2 = uj g (1 — 4& 2 ) 1 / 2 (obtained at the end of the first 
Brillouin zone k = n/ro) which is absent in the continuum limit, viz. uj 2 rj uj 2 — ujq k 2 Tq (for k <C f^ 1 ); obviously, the 
study of wave propagation in this (k < 7r/ro) region invalidates the continuum treatment employed so far in literature. 
The essential feature of discrete dynamics, to be retained here, is the (narrow) bounded TDLW ( 'phonon') frequency 
band, limited in the interval li G [(uj 2 — Aujq) 1 ' 2 , u> g ]\ note that one thus naturally obtains the stability constraint: 
uil/uj 2 g = e < 1/4 (so that wgSR Vfc 6 [0,7r/r ]). 

We needn't go into further details concerning the linear regime, since it is covered in the literature. We shall, 
instead, see what happens if the nonlinear terms are retained, in this discrete description. 
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FIG. 1: The dispersion relation of the TDL excitations: frequency u) (normalized over uj 9 ) versus wavenumber k. The value 
of uJo/ujg (~ coupling strength) increase from top to bottom. Note that upper (less steep, continuous) curve is more likely to 
occur in a real (weakly-coupled) DP crystal. 



IV. EXISTENCE OF DISCRETE BREATHERS - ANALYSIS 



We are interested in the (possibility for the) existence of multi-mode breathers, i.e. localized (discrete) excitations 
in the form: 



q n {r) = 2_j A n (m) ex.p(imujT) , 

rn— — oo 

with A n (m) = A n (—m) for reality and |A„(m)| — > as n — > ±oo, for localization. 



(8) 



A. The formalism 



Inserting Eq. (jHJ in the equation of motion J2J, one obtains a (numerable) set of algebraic equations in the form: 



P \ \ \ 

A n+ i(m) + A fl _i(m) + C m A n (m) = ^ ^ ^ A n (m 1 )A n (m 2 )A n (m 3 ) 



mi m-2 ?7i3 



- ^2 ^2 A n (m 4 ) A n (ro 5 ) , 



where the dummy indices rnj (J = 1, 2, 5) satisfy mi + ni2 + to 3 = 7714 + to 5 = to; we have defined: 

1 — m 2 Lu 2 " 



On 



2 - 



(9) 



(10) 



In order to be more precise and gain in analytical tractability (yet somewhat losing in generality), one may assume 
that the contribution of higher (for m > 2) frequency harmonics may be neglected. Eq. (JHJ) then reduces to: 



q n (t) « 2A n {l)cosujT +A n (Q). 



(11) 



Note the zeroth-harmonic (mean displacement) term, for n = 0, which is due to the cubic term (~ a, above), and 
should vanish for a = 0. The system thus becomes (for to = 0, 1): 

A n+1 (l) + A n _r(l) + CMl) = -2®A n (l)A n (0) ~ £ [^(1)^(0) + 3^(1)A„(-1)] 

a 8 

A n+1 (0) + A n _i(0) + C A n (0) = -2-A(l)A,(-l) -6^A n (0)A„(l)A n (-l), (12) 
i.e., setting A n (l) = A n (-1) = A n and A n (0) = B n , viz. q n (t) = 2A n cosujT + B n : 



i-:-r'j.l„ -2^A n B n - tL (A n Bl + 3A 3 n ] 



,/3 



-2-Ai -<S'-AiB n 



(13) 
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We see that the amplitudes A n (B n ) of the first (zeroth) harmonic terms, corresponding to the n— th site, will be 
given by the iterative solution of Eqs. (|13[) [or, of Eqs. 10, should higher harmonics m be considered]. In specific, 
one may express (|13f) as: 

a g 

a n+1 = -Cn - Cia n + 2— a n b n + — (a n b 2 + 3a n ) = /i(a„, b n , c n ,d n ) 
K+i = -d n - C b n + 2-a 2 + 6-a 2 6„ = fo{a n ,b n ,Cn,d n ) 

Cn+l = »n 

d n +l = b n , (14) 

and then iterate, for a given initial condition (ai, bi, c\, di) = [A\, B\, Ao, Bo), the map defined by (|14Jl . 

At this stage, one needs to determine whether the fixed point of the 4-dimensional map l|14f> [or of the complete 
4N-dimensional map corresponding to @, in general] is hyperbolic, and examine the dimensionality of its stable and 
unstable manifolds. It is known that the existence of discrete breathers is associated with homoclinic orbits, 

implying a saddle point at the origin. 

Let us now linearize the map Ijl4(l near the fixed point {a\, &i, c\, di) = (0,0,0,0) = O4, by setting e.g. 
(a n ,b n ,c n ,d n ) = (£1, £2) £i)n = € 5i 4 , where £j <n < 1 (j = 1, ...,4). One thus obtains the matrix relation: 



where M is the matrix: 



M 





= M 
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(15) 



(16) 



Now, it is a trivial algebraic exercise to show that the characteristic polynomial p(X) = Det(M — AI) of this matrix 
may be factorized as: 

p(A) = (A 2 + C A + 1) (A 2 + C t X + 1) = Po(A)pi(A) , 

implying the existence of 4 eigenvalues, say Ai 2,3,4) such that po(Ai,2) = Po(A3,4) = 0. One may check that the 
condition for all eigenvalues to be real and different, hence for O4 to be a saddle point, amounts to the constraint: 
|C ,i| > 2, i.e. C i [-2,2] and d £ [-2,2]. Recalling that 

Ci = (l-2e-c; 2 )/e, C = (1 - 2e)/e, (17) 

from l|10|l . one finds the (simultaneous) constraints: 1 — 4e > and (1 — lo 2 )(\ — lo 2 — 4e) > 0. One immediately sees 
that the former (i.e. e < 1/4) corresponds to the linear stability condition mentioned above, while the latter amounts 
to the requirement that the breather frequency should lie outside the 'phonon band', viz. oj 2 /u> 2 ^ [1 — 4e, 1]. 

It is straightforward to show that in case one considers the complete multi-mode map, defined by Eq. @, one 
obtains an analogous factorizable characteristic polynomial for the 4A^ x 4N matrix M, viz. p(X) = [] m Pm(A). The 
same analysis then leads to the hyperbolicity criterion: 

|C m |<2 m = 0,1,2,... 

One thus recovers, in addition to the first of the above constraint (e < 1/4), the condition: mui/uig </ (1 — 4e, l) 1 / 2 
(Vm = 0,1,2,...), which coincides with the - physically meaningful - non-breather-phonon- resonance condition found 
via different analytical methods 0, 0, 0] . We see that the breather frequency, as well as all its multiples (harmonics) 
should lie outside the allowed linear vibration frequency band, otherwise the breather may enter in resonance with the 
linear TDLW ( 'phonon') dispersion curve, resulting in its being decomposed into a superposition of linear excitations 
(and hence de- localized). 



B. Numerical analysis 

At this stage, one is left with task of finding the numerical values of A n ,B n [ci. I|13|) ] for a given homoclinic 
orbit; these may then be used as an initial condition, in order to solve the equation (|13|) numerically, by considering 
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FIG. 2: Localized discrete breather dust lattice excitations; the successive lattice site displacements are depicted at maximum 
amplitude: (a) odd-parity solution; (b) even-parity solution. 



a given number of particles N and harmonic modes m max (viz. m = 0, 1, 2, m max ). One thus obtains a given 
set of numerical values for u n (n — 1,2, ...,N), which constitute the numerical solution for the anticipated breather 
excitation. The stability of the solution thus obtained, say q n , my be checked by directly substituting with q n = q n +£,n 
(for n — —N, 0, N) into the initial equation of motion |J3}. 

This numerical scheme is now being elaborated, and the detailed results will be reported in an extended paper, in 
preparation. 



V. BREATHER CONTROL. 



The stability of a breather excitation may be controlled via external feedback, as known from one-dimensional 
discrete solid chains [3i| • The method consists in using the knowledge of a reference state (unstable breather), say 

8zn^ = z n (t), e.g. obtained via an investigation of the homoclinic orbits of the 2d map obeyed by the main Fourier 
component [9], and then perturbing the evolution equation (1) by adding a term +K[z n (t) — Sz n ] in the right-hand 
side (rhs), in order to stabilize breathers via tuning of the control parameter K. This method relies on the application 
of the continuous feedback control (c/c) formalism (see the Refs. in |l3j ). Alternatively, as argued in |l3| . a more 
efficient scheme should instead involve a term +Ld[z n (t) — 6z n ]/dt in the rhs of Eq. (1) (dissipative cfc), whence the 
damping imposed results in a higher convergence to the desired solution z n {t). Preliminary work in this direction is 
being carried out and progress will be reported later. 



VI. CONCLUSIONS - DISCUSSION 



We have investigated, from first principles, the possibility of existence of localized discrete breather-type excitations 
associated with vertical dust grain motion in a dust mono-layer, which is assumed to be one-dimensional. 

It may be noted, in concluding, that the localized structures presented in this Letter, owe their existence to the 
intrinsic lattice discreteness in combination with the nonlinearity of the plasma sheath. Both are experimentally 
tunable physical mechanisms, so our results may be investigated (and will hopefully be verified) by appropriately 
designed experiments. The experimental confirmation of their existence in dust crystals appears as a promising field, 
which may open new directions e.g. in the design of applications. 
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